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FOREWORD © 


Since changes that occur in the thermal structure of 
the oceans have a far reaching effect on various naval 
operations, any practical technique that can be developed 
to predict these changes would have considerable oper= 
ational utility. 

This report represents the results of research con- 
cerning the effect of weather on the thermal structure of 
the ocean conducted by the Division of Oceanography up to 
30 June 1952 under ONR Contract 260005. This project is 
being continued by the Hydrographic Office with periodic 
progress reports to be published at appropriate intervals. 

It is requested that activities receiving this pub- 


lication forward their comments to the Hydrographic Of- 


ye Spe 


KARL F. POEHLMANN 
Captain, USN 
Hydrographer 


fice. 


cI bal 


“ tw om v 
Bia _ pa ee 
aye 5 ei Fa 


f 
oa 


ae 


Loven suet ene ce dosti Serres see. eS Ce ae PS 


santosh we tik (oe? aie A660. doe) fates eae 


—itcm Segement Gren > 


gQoutadov ersais Yel oaee 
: i) 

pup © x ve ¥ a 

eter bi canes ine 


" 
j 


eeoo Coveeeert to. co Mines edt arhanetget 29 Mar chases 


a o Oe a: e F es - oa ae. ety ns tas Deere =f 
So Sine sees Laansedt ould ao si iYest te opts ogy aa 


oF oe ek esas re te A er ert ina Kod Breese SRO Saree 
: j 
Ste Tie = d sCUUGGS Tae it te See eS 


piawieial etabxeoniisa Ja bee blauce ao ). SQA ye: tae 
a? eo eee! ai A eS Ge ik. tan eae el ke 
ae SLi - cpeeletoy ears Oce tage BO i Bas ie sb. 


ait) ptrigesahedysd a ed 4bttihoeme a vied atawse'’ aokael 


» “ 
vr way nig 
" ia ihe 7 s 
= _ oe ¥ 
wis Ay Ft . 7 
‘ Te 
v oy > 
bas nig wae 2 
‘ ~erd 
aj t ri 
= 1 
+ 
= 
v 

i 
i I i 

i 

t i 
{ , : 
id 


ABSTRACT 


WAS jOLOOOSS OL WAL PPOJSOG LS wO WAVES wlseAwS we we wsowo— 
Vogical factors that affect materially the thermal structure 
of the ocean and to develop techniques capable of predicting it. 
This report contains: (a) a description of the factors known 
to affect the thermal structure; (b) a summary of the known 
principles concerning each of these factors; (c) proposed predic- 
tion techniques based on these principles; (d) proceedures and 
results of technique development so far undertaken by the pro- 
ject; and (e) proposals for future work by the project. 

Two types of processes were studied: (a) those that are 
caused mainly by various meteorological variables; and (b) the 
process of vertical motion caused by internal waves. Those 
processes due to the horizontal field of m>otion were disregard- 
ed because ot inadequate data. 

Based on known or predicted meteorological variables, the 
processes of radiation, evaporation and molecular conduction, 
convection, and turbulent mixing are considered. ‘The relation- 
ships civing the magnitude of the effect of these processes on 
the thermal structure are examined and typical values are pre- 
sented under mean and extreme conditions. ‘These relationships 
are fused into a proposed prediction technique. 

A sample forecast illustrated by figures and graphs is 
made by considering the above processes and using the above 
mentioned techniques. Results are promisinz enouch to warrant 
further consideration of the seneral approach. However, data 


presented indicate that an accurate and operationally useful 
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method of predicting the thermal structure cannot be obtained 
without consideration of the internal wave. 

This problem is investigated by using two models, namely, 
the two-layer system and a model assuming the density to be a 
continuous function of depth. 

Since the two-layer system assumes the wave to occur along 
the interface between fluids of different densities, the problem 


Hsp ted to the Studysos tae Oscillations of vie thermociane 


with time. It appears that the main variations in the thermocline 
depths are: (a) seasonal variations; (b) lunar or monthly vari- 
ations; (c) semidiurnal and diurnal tidal variations; (d) short 
period variations; and (e) random variations due to mixing. Data 


are presented to illustrate all these effects. 


Because previous studies have shown (c) to be the larcest 
SLMENS SEBSCG, We LS SisMelrecl Eilodwaleelidhy elie Siomile Weiaerwias iA 
Honrmulal ls) presenbed Givanes the amplitudes of these oscillations 
imeone area as a function of thie thicknesses of and density dis= 
continuity between the layers by assuming a two-layer system on 
a rotating spherical earth. In the few cases in which the data 
were adequate for a forecast, the relationship predicted the 
amplitudes with a fair degree of accuracy. 

When the continuous density model is considered, a more 
Conplrcarecd pacLuUre) of ithe winternak oscilla tvons resuilius am 
which the amplitudes and phases vary continuously with depth. 
Methods are developed for predicting these internal oscillations 


Zac 2 Seimgolls jucscliGipiCia as jomsseaigeel. Wale OMOsiciayetel yall joresicliS cel 


amplitudes show good agreement at the greater depths, but near 
the surface discrepancies occur. Possible reasons for these dis- 
ecrepancies are examined;it is proposed that deficiencies of the 
available data comprise the largest single factor. 

This report was prepared by J. J. Schule, Jr., with the 


assistance of L. S. Simpson and A. Shapiro. 
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Temperature changes over 3 ‘days due to loss of 
heat from oceans to atmosphere. 


Change in thermal structure due to turbulent 
mixing. 


Pyrheliometer record, 1600-13002, 22 April 1951. 
BT trace, 0500Z, 20 April 1951. 
Transparency station, 23 April 1951. 


Energy spectrum of radiation from sun and sky 
at the surface. 


Change in temperature distribution due to 
radiation, 20-27 April 1951. 


Change in tempereture distribution due to 
heat exchange with atmosphere and due to 
convection, 20-27 April, 1951. 


Change in temperature distribution due to 
mixing, 20-27 April 1951. 


BT trace, 2000Z, 27 April 1951. 


Generation of internal waves by an undersea 
ridge. 


Mean monthly mixed layer depths for weather 
ship COCA. 


Mean monthly mixed layer depths for weather 
ship DELTA. 


Mean monthly mixed layer depths for weather 
ship ECHO. 


Mean monthly mixed layer depths for Pacific 
weather ship NECTAR. 


Mean monthly mixed layer depths for Pacific 
weather ship GOLF. 


Mean monthly mixed layer depths for Pacific 
weather ship PAPA. 
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I. Purpose 
To investigate the effect of weather on the thermal 


structure of the ocean, 


II. Procedure 
A. Introduction 

This project is aimed primarily at the eventual 
development of practical and operationally useful techniques 
for the prediction of the thermal structure of the ocean. 
Unfortunately, this is a complicated and difficult problem, 
for the thermal structure is affected by many processes, 
some of profound effect and all seeking to change its nature 
simultaneously. Separating these processes for study and 
estimatin= quantitatively their relative effects is therefore 
the first field of interest for this project. Added'to’ the 
difficulties inherent in the problem itself is the fact that, 
in a science where controlled experimentation is impossible 
and data collection expensive, very little data adequate for 
an investigation of this type are available. 

We may consider that the processes that affect 
the thermal structure may be divided into three classes: 
the thermodynamic, which involve a change in the heat content 
of the water column; the static, which involve temperature 
changes in the column due either to internal or external forces 
without change in the heat content; and the dynamic, which are the 
result of the three-dimensional field of motion. Under each 
we may list several processes: 

1. Thermodynamic processes 
ae Radiation - incident, reflected, and back radiation. 


b. Exchange with atmosphere by molecular conduction. 


ec. lxchange with atmosphere by evaporation. 
2. Mixing processes 

a Convection 

be. Turbulent mixing 
3. Dynamic processes 

a. Horizontal advection 

b. Vertical advection due to internal waves 

The parameters that affect the thermodynamic and mixing 
processes are mainly meteorological, such as air temperature 
and humidity, vapor pressure, wind, sea and swell, etc. It was 
therefore decided to attempt, by consideration off known principles, 
to evaluate the effects of these processes first, at the same time 
Peeemp tine to minimize the dynamic effects by testing these 
relationships with data collected in areas where large scale 
moptons Go Not occur. lh this atitempy were successnull, the 
dynamic effects then could be considered. 

The following discussion should be considered more as an 
exposition of possible methods and techniques that might be 
menitzod in attacking this problem than as a body of proven 
scientific results leading to a solution of the problem. 
Although several of tne procedures tested yielded promising 
results, the few months that the project has been under way 
and the lack of available data have precluded the opportunity 
of testing a sufficient number of cases to demonstrate their 


validity conclusively: 


B. Direct investigation of meteorological parameters 
1. Introduction 
We might divide the processes that depend on 
meteorological parameters and affect the thermal structure 
of the oceans into two classes: (1) the thermodynamic processes, 
which have to do with the actual addition or subtraction of 
heat from the ocean, and (2) the mixing processes, which cause 
changes in the thermal structure due to some internal or 
external force on the water column without changing the actual 
heat content of the column. The dynamic processes, which are 
those due to the three-dimensional distribution of motion in 
the oceans and the forces involved,have been neglected in this 
portion of the report. 
2. Thermodynamic processes: Radiation 
ae incoming radiation 

The most complete treatment of the energy exchange 
between the atmosphere and the ocean has been given by Jacobs 
(1942). Although this paper approached the problem from a 
climatological viewpoint, the relationships are of such a nature 
that, by knowing the values of the constants for a given area 
and being able to predict the values of the meteorological 
parameters involved, the total amount of heat gained or lost by 
the ocean during a given period may be determined. It is then 
necessary to translate this amount of heat into temperature 
changes; here the problem is complicated by the fact that we are 
interested not in the surface temperature change alone but in the 


change at every depth, 


It is necessary to consider first the total amount of 
heat that reaches a unit area of the sea surface in a given time 
and is thus available to the water column for absorption at 
various depths. Mosby (1936) has shown that, with clear skies, 
the total incoming radiation may be accurately approximated 


by the simple linear equation 


Q¢ kz (a) 


Here Qo is the total incoming radiation expressed in g.cal. 

emr@ minzl, x may be considered as a parameter whose value 

varies with the water vapor in the atmosphere, but its variation 
is in such a narrow range that, for the purposes of this study, 

it may be considered as a constant and equal to 0.025. ais the 
mean solar altitude, in degrees, during the period over which the 
heat is calculated. 

This expression must be modified for cloudiness. This 
problem has been treated empirically by Angstrom (1922) and 
Kimball (1928). ‘Their findings indicate that the values for 
the ration of overcast sky radiation to clear sky radiation 
vary between 0.22 and 0.30, the most likely mean value being 


0.29. Applying this ratio, we have 


C 
Q,=.0252[0.29+0.7I(I- 55)] (2) 


Here Q; is the total amount of incoming solar energy reaching 
the sea surface in gs. cal. emz2 minz+ and © is the percentage 
of sky covered by clouds. 

It is now necessary to apply a correction factor to equation 
(2) to account for the percentage of the incoming radiation 
that is reflected back to space. Schmidt (1950) has determined 
this percentage to vary from 3.3 at the Equator to 8.0 at the 
poles. Designating this percentage as r and applying it to 
equation (2), we have a final expression for the total amount 


of incoming radiation, namely, 


‘4 C 
Q,,7.025z [0.29+0.71(I- Ra) (I-n) (3) 


Qap is still expressed in g. Gals ema wile if it 2s depieegec 
to calculate the total incoming radiation for some period of 
time, it is only necessary to multiply Q,, by the period t 
expressed in minutes. 
195 Back radiation 

Not all of the heat energy defined in equation 
(3) i8 actually absorbed by the water mass.. Some of this heat 
will be radiated back into space. -Angstrom (1920) has defined a 
quantity, eff. Q,), which depends on air peuroeen EUS and humidity 


and gives the back radiation with a clear sky. These values are 


tabulated in Angstrom (1920) and are presented in nomogram form, 
slightly modified for a sea surface, in The Oceans (192).  “eeiss 
value also must be adjusted for the amount of cloudiness present 


and for internal reflection. When the most generally accepted 


values for these corrections are applied, the following equation 


results: 


Q,=0.94[eff Q,(I-0.0083 C)] (4) 


In this expression Qp is the average back radiation in 
a. cal. emz< minz+ and G is the mean cloudiness in percent. 
Just as in the incominge’radiation, for a finite period of time 
it is necessary to multiplty Q, by the period t expressed in 
minutes. 
ec. Change in temperature with depth due to absorption 
of radiation 
Subtracting Q, from Q,, will yield the amount 
of heat available for heating the oceans and for exchange with 
the atmosphere. To proceed with the investigation it was 
necessary to make two assumptions: (1) for sufficiently short 
period of time the processes of heating the water and of exchange 
of heat between ocean and atmosphere may be considered mutually 
exclusive processes and may be computed as such, and (2) each 
layer radiates in proportion to the radiation received. 
Neglecting for the present the process of exchange with 
the atmosphere, therefore, we may write an expression for the 


temperature change at any depth, namely, 
GS) 
AT a Cel 


where 
Al(zZ)= temperature change at depth z, 


Gua eS ecectlm ce heat or Wauer ab) cons tanh pressure 
(approximately 0.93) 


7 


There remains the problem of evaluating f(z). This may be 
done in one of two ways. If the extinction coefficients for 


the area are known, we may employ the relationship 


-kZ 
(Z)= 
f(Z)=e (6) 


Here k is defined as the extinction coefficient. For practical 
use it has become the convention in oceanography to consider k 
a constant over a meter's thickness; therefore, k is generally 
given per meter. For purposes of computation it is also more 
convenient to utilize equation (6) over succeeding intervals of 
one meter or more; if one meter is used, 2 = 1, if two meters 


are used, Z ms 2, etc. 


Now, if we start at the surface with (Qab -&p) and designate 
this quantity as Q,;, the amount of incoming heat left at the 


bottom of the first meter is 


my -k 
Pe (7) 


where kj is the extinction coefficient in the first meter. 


The amount of heat absorbed in the first meter is 


A Q,, = QF Or e 


By continuing the process, the amount of heat remaining at the 


pottom of the second meter is obtained. 


-k, -k (Kk, +Kp) 
Q,,7Q,e 'e 7=Qre | ? & 


where kp is the extinction coefficient between the first and 
second meters. The amount absorbed in the second meter is 


v Ki (Ky tk) 
AQ,,=Q,[e ena 2] (10) 


In general, the amount of heat absorbed in the nth layer is 


n-t n 
“2k ~ Zk; 
AQ,,=Qr le Fo =a | (11) 


where k=O, 


Now, dividing the values of AQ, for the various layers 
by the thickness of the layers in centimeters and by cp 
will yield the temperature change in that layer in °C. 

A simpler method of computing temperature changes at the 
different depths sometimes may be employed. When transparency 
stations are occupied at sea, the data may be collected in 
terms of the percentage of incoming radiation at the surface 
that is left at any given depth. If Q, is known, therefore, 
the AQ, for any layer may be computed directly from these data 


rather than from the extinction coefficients. 


d. Typical values of temperature changes at 
different depths due to the absorption of 
radiation 

The above considerations have been used to compute 
the data presented in table I, which gives an indication 
of the magnitude of the temperature changes that might be 
expected at different depths and at various latitudes resulting 
from the absorption of radiation alone. These data were 
computed for a typical February and August day by using a rea- 
sonable mean altitude of the sun in each case. The percentage 
of cloudiness was assumed to be 50, and r was allowed to vary 
linearly with latitude within the limits described in section 
2a above, the extinction coefficients were taken from a table 
of typical values for clear oceanic water presented on page 


107 of The Oceans (192) 


Table 7 


Total Heat 

Added (g.cals| Temperature Change °F 
late| cme day ee ce on) (wert T00=103 

eee 318 | 2.29|2.26].10 | .10| .02 007 
30 orn] 52 | heer 16188) 2:¢o\Woueh’ a1 soul eee 
70°N| -1.45] 314 = [2.03 Se estore). ) 


Ly 
e 
ine) 

i 
e 
QD 
(e) 
(oy 


Typical daily temperature changes (°F)due to 
absorption of radiation at various latitudes 
ane depth intervals. 
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3. Thermodynamic processes: Exchange of heat between 
ocean and atmosphere 
As has been discussed above, all the absorbed heat 
does not go into heating the oceans; some of the heat is lost 
from the water, because of exchange with the atmosphere. This 
exchange is composed of two processes, evaporation and molec= 
ular conduction. According to Jacobs (19:2), both these 


effects may be incorporated into the following single relation= 
ship: 
ToT, 
Q,145.4Wiley-e,)[-01( BE 2} | 


ew &g (12) 


where 


Q,= amount of heat gained or lost by the sea surface 
einweny Caluaren ae 


W =wind velocity in knots, 
t 


time in days, 


€y= saturation vapor pressure over the water surface 
alias ale We 


@,= vapor pressure of the air (in. Hg.), 

Ty= temperature of the water in °F., 

T,= temperature of the air in ory 

An indication of the orders of magnitude of Qe under 

typical conditions may be obtained from table II. These 
data were computed under the assumption that the air temperature 
was 60°F, and the relative humidity 75%. Values for Q, are 
presented for various wind speeds ana for various sea-air 
temperature differences by assuming that the process of exchange 


went on for one day only. it will be noted that with high 


ll 


winds and large sea-air tue datierences, such as smiehir 
be Found don tanee a strong polar front, Q, is extremely large 
and misht easily be greater in magnitude than Q@p. When this 
occurs the net flux of heat per day is negative with respect 
to the ocean; therefore, to arrive at the new thermal structure 
it would be necessary to consider the process of convectional 
cooling to relatively great depths. If, however, the flux 
or heat is positive with respect to the oceans, the conveetvon 
resulting from exchange with the atmosphere will take place 
within a relatively shallow layer because of marked surface 
stability. This is the reason why winter temperature traces 
are essentially isothermal to great depths whereas a shallow 
thermocline exists in summer. Superimposed upon this effect 
is the effect of turbulent mixing, which will further modify 
the thermal structure. The processes of convection and 
turbulent mixing will be discussed in the following paragraphs. 
lh. Mixing processes; Convection 

So far we have developed the following constituents of 
a predicted thermal structure: (1) an amended temperature 
distribution, T,(z), which would result from a process of 
radiation alone (discussed in 2 above), and (2) an amount of 
heat, @g, which is to be lost to the water because of exchange 
with the atmosphere during the forecast period (discussed in 
3 above). 

Our problem now is to determine the modifications that 

WallLdl exer@iwhe aha WAG) due to the loss of heat Qe at the surface. 


When one considers this problem, it is necessary to consider 
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the effect of salinity on this new thermal structure. Prac- 
tically, however, the only available data may consist of BT's 
with no accompanying salinity data. Therefore, it would be 
wise to consider possible methods for either eventuality. 


ae If salinity data are not available in this case, 


we may write 


h 
Q,=c,f[T, (2-T(h) Jaz 


PJ, (13) 
where T(h) is the temperature that exists from the surface. 
to some depth h which must be chosen to satisfy the integral. 
Since the value of this integral, Q,, is known, if Tp(z) is 
known explicitly, then T(h) may be easily calculated. However, 
since we do not always know the function explicitly, it is 
much more convenient to evaluate this integral graphically. 
nanee the units of Eye UGE | SlatvieTa al fe6) eleMe 3 om>*, the area, 
[T, (2) ~7(h)]dz must have the units °C cm. Since BT data is 
given in °F, and feet, this area may be numerically evaluated 
by a simple conversion factor 
h 
Q,=16.95¢, f [J zy-T(h)] dz 
e (1) 
if it becomes convenient to’ use meters’ for depth instead of 


Heet in plotting the curves, the conversion will then be 


h 
Q,=55.6 ea La (z)-T (h)| dz 
(o) 


An indication of how some of the typical values of G,& 
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as given in table II may be applied to different types of 
Leniperauire Gistripattons is presented in figure 1, Here curve 
A might represent Typ(z) (typical temperature change) for winter, 
and curve B for summer. To show the effect more clearly, Qe 

was taken over a three-day period. In this figure the effect 
of the initial temperature distribution is noticeable. With 


a typical shallow summer thermocline, as is shown in curve B, 


Table II 


WIND SPEED (KNOTS ) 


i 


2 


Amount of heat (gs. cal. cm. day ) wha 16 los co 
the ocean due to exchange with the atmosphere for 

various wind speeds and sea-air temperature differ- 
SiNGSS - 2 


the change in temperature of the mixed layer will be much 
sreater than with a deep thermocline for an acne Loss) ol 
Heat. im addition, it will be noticed IAG WAG) Smee eI Ot 
more mPence the discontinuity at the thermocline the less 


will be the change in the depth of the thermocline due to 


convectional cooling. 


Oo, Le Sedliianiy Glensel cues eiyetlelolie sal adadition to 
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temperature data, more accurate methods of predicting the thermal 
structure resulting from a given loss of heat, Qe, at the 
surface are possible. Such methods already have been devised 
by Zubov (1933) and Defant (1949). Detailed descriptions and 
examples are given in the references and, since the methods 
were not employed in any of the examples used in this report, 
will not be repeated here. Essentially, the technique involves 
dividing a column of water into several small layers and, 
by successively mixing deeper and deeper layers, computing a 
new distribution of temperature, salinity, and density, as well 
as the amount of heat necessarily lost to the ocean to bring 
about this convection. In this way, it is possible to stop 
the process at the layer where this amount of heat loss is 
equal to the given amount Qe. 
5. Mixing processes: Turbulent mixing 

We now have arrived at a predicted thermal structure 

which includes the effects of radiation, exchange of heat with 


the atmosphere, and convection. Let this temperature distri- 


bution be called Tg(z). Depending on the amount of convection, 
this temperature distribution may or may not include a mixed 
layer, which may be called Dj. The remaining problem is to 
determine whether there will be enough turbulence generated by 
the stress of wind or sea to cause the water to be mixed to 
a depth greater than Do. 

No adequate practical techniques for attacking this pro- 
blem exist in the literature. However, a possible approach 


to the problem was presented in an earlier progress report 
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of this project (inc. Hydro ltr code 3-JH/mlb serial 336). 
This approach is partly deductive, partly empirical in nature, 
and assumes that the depth to which mixing will occur is a 
function of three parameters: 
a. The amount by which the predicted wind speed will 
exceed the critical wind speed Ve as demonstrated by Munk (197). 
b. The wave number, 8 ,associated with the predicted 


sea state (8=*7) j 


ec. The initial stability age » designated by E. 
Under these assumptions the following relationship was 
derived: 
pi | - 2 
(1+ -)(I+ 2D) ? oy 
where 


D = depth of thermocline due to mixing, 
Dj= original depth of thermocline, 


A,K = constants of proportionality which must be determined 
by data, 


density in mixed and unmixed layers, respectively, 


PoP 
t =time (units depend on units chosen for V and Vg). 
Further work must be done to determine whether this 
relationship can be developed as a practical prediction tool. 
For this reason, prediction examples have been confined where 
possible to those cases where the mean value of V was less than 
the 13 knots given as the value of V, by Munk. 


Once the mixed layer due to turbulence D has been arrived 


uy 


at, it is necessary to amend T, (z) to inelude this deeper 


thermocline. 


The temperature that the water will have through 


this mixed layer may be represented as 


1c? 
ips Bh T, (2) dz 


It is once again more convenient 
by graphical methods rather than 


T,(2) precisely. Figure 2 is an 


TEMPERATURE 


ADJUSTMENT TO ACCOUNT 
FOR TRANSITION ZONE 


DEPTH 


PREDICTED SOUNDING PRIOR TO MIXING 
PREDICTED SOUNDING AFTER MIXING 


WWUE: AREA OF COOLING QUE TO MIXING 
K WN AREA OF HEATING DUE TO MIXING 


FIGURE 2. CHANGE IN THERMAL STRUCTURE DUE TO TURBULENT 
MIXING. 


(lm 


to evaluate this integral 
determine the function 
illustration of how this might 
be done. Since no heat is add- 
ed or subtracted from the water 
column aS was the case in Sec. 
| above, the area representing 
the heated portion should equal 
the area representing the cool- 
ing portion. A fair approxima- 
tion may be made simply by 
counting squares. It should be 
mentioned that in this process 


the assumption must be made 


that the discontinuity between 


the new mixed layer D and the former sounding De) is of zero 


order. 


Adjustment must be made for this, since it is known that 


the thermocline is not an ideal discontinuity but that a tran- 


Sition zone exists. 


This adjustment would have to be made in 


accordance with continuity of the past soundings and experience, 


which automatically brings in a factortiof "art," as opposed Ga 
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"science," 


similar to that existing in modern weather fore- 
Castine in hielne ja they tinal forceacjbed sounding is! shown 
Piuite sold Vane, 
6. Example of forecast 
ae Introduction 
The methods for estimating the effect of the meteor- 
ological parameters on the thermal structure, as outlined in 
Secs 2 to 5 above, were tested by attempting a prediction 
utilizing actual data. For this purpose there was available 
a collection of observations taken by the U. S. Navy Hydrographic 
Office on cruise AMOS VIII in April 1951. On this cruise, some 
of the latest instrumentation developed by the Oceanic Develop- 
ment Branch, Division of Oceanography, were employed. It was 
Secuced tO Ubalize Lor the forceast that portion eh the ldata 
gathered in the Sargasso Sea area, since here the data were 
more complete and the mean wind speeds were generally 12 knots 
or below, thus allowing the assumption that the turbulent effect 
was small. These data include: (a) hourly BT's; (b) hourly 
observations of meteorological variables; (c) hourly observat- 
ions of sea and swell; (d) transparency stations using photocell 
equipment every few days; and (e) continuous pyrheliometer 
recording of incoming and outgoing radiation. 
The type of radiation data available is illustrated 
in figure 3. It will be noted that these data were taken 
by three separate pyrheliometers, each recording directly 
BiemCnne Of radiatvon an fecal. emz> minzl, Pyrheliometers 


1 and 2 were mounted on either side of the ship to measure 
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back and reflected radiation from the sea surface; to arrive 

at the total effect the reading from the two instruments must 

be added. Pyrheliometer 3 measured the total incoming radiation 
from both sun and sky. 

It was decided to attempt to forecast the changes in the 
characteristics of the thermal structure from the morning of 
POrApral to the vevening of 27 April, A representative BI, 

taken at 0500Z 20 April, was chosen 


TEMPERATURE (°F) 
poe com c So OMS CIO as Lewin al tompera buce (distr. ou bom. 


This BT trace is reproduced in figure 


l and plotted in meters for the sake of 


convenience. 


It also was decided that the pre- 


diction of the meteorological variables 


themselves was not within the scope of 


DEPTH (METERS) 


this problem, Therefore, the meteor- 


ological values as actually observed 


on the cruise were used in making the 


oceanographic forecast. Table III 


presents this meteorological data for 


FIGURE 4. BT TRACE, 05002, 20 APRIL 195. Gach day, together with the mean daily 
eiurcude Of the sum and the number of hours of sunlight. The 
values for the sun's altitude were extracted from tables in 
EMemebaaiaimanac lor WOSM land i sO4 Pubes No. 208, by assuming a 
mean daily hour angle of 45° during the daylight hours. 

b. Procedure 


(1) Radiation (see Sec. 2 above) 
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The initial step was to calculate for each day the 
average value of Qap according to equation (3). These values 
for the total incoming radiation are given in the first column 
of table IV. In like manner, values for Qp may be computed 
by using meteorological values given in table III and equation 
(4). These computed values for the back radiation appear in 


the second column of table IV. Both these sets of values are 
2 aL 


Seems ICAL em. aiming ys 
Table IV 
Computed Observed 
Date (g.cal.cmr2minzt) (g.cal,omzemin,~1) 
Qab ea Qb 23.0 “nef 
fee 605 61 lowobagy 
[at see | a init 

Afeulis ee, 
551 6 LOO 
»633 
667 

26 Apr 910 . 896 

2 i Apr - 760 ~Li2 5 31.0 


*Inaccurate pyrheliometer due to shock. 
Computed and measured rates of total incoming 
and back radiation (g. cal. em7< min.-l) for 
20-27 April 1951. 
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The first problem was to determine how closely the computed 
values for incoming radiation represented actual conditions dur- 
ing the forecast period. To do this pyrheliometer records were 
neeces For each day during the daylight hours the average rate 
of incoming radiation was computed from the curve of pyrhelio-~ 
meter record number 3, by using finite differences at 15-minute 
intervals. These observed values are given in the third column 
of table IV. A similar process was carried out with curves of 
pyrheliometer record numbers 1 and 2. Adding the two gave the 
average daily rate of reflected radiation. These values for 
Qref are shown in the fourth column of table IV. The similari- 
ties between the computed and observed rates of incoming radi- 
ation, -as shown in the first and third columns of ‘table IV, are 
remarkable 

The next step in the computation was to get the total 
incoming and back radiation for the eight-day period. Adding 
the values for Bais in‘the first column of table IV and multi- 
plying by 13 hours or 780 minutes give a total ine oman rad- 
iation of 607 g. cal. emr2, Adding the values of Q, in the 
second column of table IV and multiplying by 2), hours or 110 
minutes give the total Gin cow oe "cals. eames Subtracting Qy 
from Qab gives a total of 3,339 g. cal. ems absorbed by the 
water column during the period. 

It is now necessary to convert this amount of absorbed 
heat to temperature changes with depth. A representative 
transparency station, occupied on 23 April, was chosen. Vests 
data were taken by the use of colored filters, thus giving a 


separate curve for each spectral band, Figure 5S presents the 
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SURFACE LIGHT (%) > 
ion Cava for the red, yellow, 


the spectra, For each of 


these bands the graph shows 


the percentage, on a log- 


DEPTH (METERS) 


arithmic scale, of the in-= 


coming radiation that is 


7 green, and blue portions of 


FIGURE 5. TRANSPARENCY STATION. 23 APRIL 1951. left at any depth. The 
problem is to determine how much of the radiation is concentrated 
in each band. Figure 6 is an adaptation of a graph presented on 
page 105 of The Oceans (192) showing the distribution of the 


energy spectrum at the surface of the radiation from sun and sky, 


a 
fe) 
[e) 


plotted as a function of wave length. 


TS 
fo} 
(eo) 


Taking 0.64 micron as the lower end 
of the red band and measuring the area 


under the curve, one arrives at a fig- 


ENERGY IN ARBITRARY UNITS 


APEC Neos Ti micas : ure of 64.3% for the red and 35.7% for 

ee ae the remaining bands. It was decided 
to use the curve for the green band as representative of the re- 
maining 35.7%. 

By starting with the total absorbed heat of 3,339 g. cal. 
emz© it now can be determined that Coun? @. call. emse were 
absorbed from the red band and 1,192 g. cal. cms2 were absorbed 
from all other bands. By using these figures in conjunction 
with the appropriate curve in figure 5, table V was constructed, 
showing the amounts absorbed from each band in each layer. The 


totad number of g. cal. cmz2 absorbed in each layer for all bands 
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is shown in the extreme right hand column of table V. 

To convert these totals of the amount of heat absorbed 
in each layer into temperatures, it is only necessary to divide 
by the number of centimeters in each layer and by Cp (.93). 
These temperatures will then be in °C; for use with BT data 
it is necessary to convert them to °F. Table VI gives the 
results of this computation with °C. in the left-hand column 


and °F, in the right-hand column. 


Tab ey a 


Layer 
| (meters ) 


Temperature changes in the various 
layers due to the absorption of rad- 
iation during the period 20-27 April 
1951 


The final step is to apply the temperature changes given 


in table VI to the original trace as shown in figure 4. ‘The 
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FIGURE 7. CHANGE IN TEMPERATURE DISTRIBUTION DUE YO RADIATION, 20-27 APRIL 1951. 


resulting temperature distribution is T,(z), and is shown in 
figure 7. 
(2) Heat exchange with atmosphere (see Sec. 3 above) 

It is now necessary to compute the amount of heat 
gained or lost by the ocean from exchange with the atmosphere. 
This may be done easily for each day of the forecast period 
by the use of equation (12) and the meteorological data present- 
sil ale) ieee wye; IWITIhe) Ulavey Tete Tope Q, for each day are presented 
in table VII. Addition of these daily values of Q, gives a 
total of 2,322 2. cal. om7*, which is the total heat lost to 


the ocean by exchange with the atmosphere. 
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Table VII 


Qe : 
(g.cal.cm-2day~1) 


219.33 


207.06 


Daily heat totals (ga cal. cms2) dost 
to the ocean, due to exchange with the 
ieee ae for the period 20-27 April 
1951. 


(3) Convection (see Sec. above) 


T(z) now must be adjusted to take care of the total 


al 
heat, Q,, lost to the ocean as computed above. Since figure 
7 is plotted in °F. and meters, it is necessary to employ 


equation (15) in this computation. Substituting for Q, the 


total given above and solving for the integral, we obtain 


uIZoeeie 


h 
Beata dA (z)-T(h)|dz=44.92 °F meters (13) 


fo) 
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Here 4.92 is the number af squares, bounded by one meter and 
one °F., which must be subtracted graphically from T,(2) to 
arrive at the predicted temperature distribution. This is done 


in figure 6, and the resulting distribution is T,(2). 
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FIGURE 8. CHANGE iN TEMPERATURE DISTRIBUTION DUE TO. HEAT EXCHANGE WITH ATMOSPHERE 
AND DUE TO CONVECTION. 20-27 APRIL 195! 


(4) Turbulent mixing (see Sec. 5 above) 

The character of T,(z) was so unusual that it was 
considered probable from past experience that the wind and 
sea would be of enough force to mix the shallow layer in 
which two abrupt discontinuities are shown. Some of this 
abruptness undoubtedly is due to the discontinuous manner in 
which the temperature changes due to radiation must be plotted. 
z) was adjusted by graphic means according to the 
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FIGURE 9. CHANGE IN TEMPERATURE DISTRI- 
BUTION DUE TO MIXING, 20-27 APRIL 195}. 
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FIGURE 10 BT TRACE, 2000Z, 27 APRIL 195 


principles of equation (17) and figure 
2. Figure 9 shows the results of this 
adjustment; the final forecast sound- 
ing for the evening of 27 April is 
shown by the solid line. Figure 10 
shows a BT observation which may be 
used as a verification of this fore- 
cast, having been taken at 2000Z on 

27 April. Except for a slight dis- 
crepancy in surface temperature, the 
agreement between the predicted and 
observed temperature distribution must 
be considered good. 

If the depth of the thermocline 
as predicted above would remain con-= 
stant for a reasonable period of time, 
it could be coneluded that a possible 
method of predicting the thermal struc- 
cure had been developed. However, 
another factor remains which requires 
investigation. It is noted that the 
thermocline oscillates in a periodic 
manner which in no way can be attri- 
buted to the effect of meteorological 
variables, It was considered that 
this problem must be considered before 


adequate means of predicting the thermal 
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structure could be developed. These oscillations, the so- 
called internal wave, are the subject of the following section 
of the report. 

G. Investigation of periodic variations due to internal 
waves 

1. Introduetion 

The problem of periodic oscillations of the water 

masses in the oceans has been treated by various authors in the 
past, and more and more attention has been given this phenomenon 
in recent years with the increasing realization of the important 
role these waves play in the determination of the thermal struc- 
ture. The problem of isolating these waves from other processes 
operating in the ocean is not a simple one; one of the facts 
that has become evident to all investigators of internal waves 
is that these waves complicate the actual motion greatly and 
cause extremely varied and intricate distributions of currents 
and vertical motions. 

In general the problem has been approached in two ways. 
The first is the simpler linear theory first treated by Stokes 
(Lamb, 1932) in which it is assumed that the wave occurs along 
the discontinuity surface between two layers, each homogeneous 
and scach of (ditterens density. The more womnplexiprebtemmou 
progressive internal waves in a heterogeneous medium in which 
density is a continuous function of depth has been attacked 
with considerable success by Fjeldstad (1933) and others. In 
oceanography both models have their applications, and therefore 


Doth have been -considered in this investigation. 
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Before proceeding with a discussion of how the two separate 
models were approached, it mizht be well to mention that the 
general conclusion reached by most investigators is that the 
sreatest single factor in determining the nature of the thermal 
structure at any one time:is the internal wave of tidal peelea. 
Consequently most of the investigations have centered around 
waves of semidiurnal and diurnal periods. 

2. The two-layer system 

This problem resolves itself mainly into a consider- 
ation of the depth and the amplitudes of the oscillations of the 
thermocline, because of the fact that the maximum amplitude of 
the wave occurs at the interface between two layers of differ- 


ent densities. It alsos can be shown, from the boundary conditions 


and the fact that density is constant through each layer, that the 
amplitude must decrease linearly from the interface and approach 
zero at both the free surface and the bottom. 

Although the problem of the internal wave was first 
attacked by Stokes (Lamb, 1932), who derived a good many of 
the relationships, the subject did not excite much practical 
interest until the work of Ekman (1904). Ekman used the concept 
of the internal wave to explain the peculiar phenomenon of 
"dead water" experienced by slow moving sailing vessels and 
steamships of that period. He showed that when "dead water" 
Was experienced, there existed a layer of light, fresh water 
overlying a heavier, more saline water below. The enersy 
which the ship usually used to overcome the resistance of the 


water itself went instead to the generation and maintenance of 
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an internal wave along this boundary, with a resultant hind- 
rance of the forward motion of the ship. 

Vertical displacements of water masses may be easily 
observed by repeating, at short time intervals, oceanographic 
duganpaticns at one place. The first investigation of this 
type was made by Helland-Hansen and Nansen (1909) on the 
Michael Sars expedition. During the cruise repeated serial 
observations were made by two ships and a harmonic analysis 
was made of the data. 

If such observations are repeated for a long enough time, 
it is possible to find by harmonic analysis the period of the 
waves. Although short periodic internal waves such as were 
described by Ekman are known to exist and are actually observed, 
the Michael Sars data yielded the information that the predom- 
inant periods of the oscillations were of semidiurnal and di- 
urnal nature. 

This conclusion has been borne out by the work of other 
authors. One of the most complete analyses of such data was 
made by Seiwell (1939, 1942). The results of this investigation 
showed that, after the oscillations of 12= and 2li-hour periods 
were eliminated, the residue exhibited a normal or random dis- 
tribution. Seiwell also demonstrated that, not only did the 
diurnal and semidiurnal waves predominate in period, but they 
also made the greatest contribution to the amplitude of the 
WAVES. 

Although these investigations have demonstrated that the 


vertical oscillations in the oceans are dominated by waves of 
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12- and 2h-hour periods, the difficulty has been to show that 
these waves have their origin in the tidal forces. Defant 
(1932) and others have shown that the tidal force is not suffi- 
cient to account for the observed amplitudes of these waves 
directly. He suggests that these oscillations may not be due 
to true internal wave motion at all, but rather are the result 
of periodic variations with time of the horizontal tidal currents, 
which cause periodic variations of the inclinations of an exist- 
ing surface of discontinuity. 

Another objection to the thesis that these waves are of a 
tidal nature has been the fact that when the velocity of propa- 
gation of these waves was computed, it was found to be much less 
than the velocity of the surface progressive tidal wave. There- 
fore, it was concluded that no resonance between the tidal forces 
and the responding medium could occur. Since the vertical component 
of the tidal force is not sufficient to causes the wave directly, and 
since no resonance could take place, it seemed to some highly unlike- 
ly that these waves had their origin in the tidal forces. 

This problem has been attacked by Haurwitz (1950), who ate 
tempted to show that resonance could occur if one considered the 


° 


effect of the earth's rotation. When this factor is considered, 


one arrives at free wave periods which are much closer in mag-= 
nitude to the semidiurnal and diurnal periods. Haurwitz also 


suggests that, if the procedure were altered to include consid- 


eration of a rotating 


) 


spherical earth rather than a rotating 
dise, much closer agreement might be obtained. This is still to 


be shown; the present investigators have made an attempt, but 
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the results in their unsimplified form are too unwieldy to be 


of much practical value, 


(See appendix A.) 


Another interesting treatment has been presented by Zeilon 


(ARSE) RS 


In this paper he stresses the importance of the bottom 


configuration in generating and maintaining oscillations in the 


Oceane 


encounters an elevation in the bottom topography. 


He discusses a model in which a uniform steady current 


He shows that 


under these conditions a wave will be generated downstream from 


the elevation and, from Continuity considerations, will also ex-= 


FIGURE || GENERATION OF INTERNAL 
WAVES BY AN UNDERSEA 
RIDGE. 

Drasran  OGeseries. Ol Sax 

photographs taken at in- 

tervals of about one-half 


second, showing the gen= 
eration of waves in the 
boundary above the ridge 


and their propagation on 
both jsiides) jot at.) ‘Touche 
left of the ridge is vis- 
ible a small pendulum in- 
dicator, moving with the 
water and showing the dir-= 
ection and intensity of 
the current. 


tend upstream. The greatest in-= 
ternal wave amplitude will be in 
the vicinity of the elevation. 
Pigure 11, taken from Zeilon's 
work, illustrates this effect. 
Zeilon's work is mentioned to sug= 
gest that the bottom contours play 
some part in the formation of in= 
ternal waves. It is therefore ad- 
visable in any study of this phen= 
omenon either to take this effect 
into account or to eliminate it, 

It was decided in this study to 
utilize data taken in one locality 
in order to eliminate bottom topog-= 
raphy effects. 


In discussing the two-layer 


system, mention should be made of a 
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new approach formulated by Davies(1951). This paper has not as 
yet been studied thoroughly, but in the future the practical 
applications of this work to the present problem will be 
evaluated. 

Since the basic assumptions made in the simplified two- 
layer model predicate the wave to be traveling along the inter- 
face, a study of the internal wave using these principles must 
moimiv be a study of the oscillations of the thermocline. One 
may begin by defining the processes which, at any given time or 
for any period of time, will govern the depth of the mixed layer. 
IMA SE eheeye 

a. seasonal variations - Probably more is known about this 
parameter than any other because it is strongly reflected in the 
mean monthly values of the thermocline, which mean values have 
been published. Tillustrations 
ef this effect, taken from these publications, are given in fig- 


ures 12, 13, and 14. Since most of the time-series BT data are 
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FIGURE 12. MEAN MONTHLY MIXED LAYER DEPTHS FOR WEATHER SHIP COCA, 52.7°N.,35.7°W. 
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FIGURE 13. MEAN MONTHLY MIXED LAYER DEPTHS FOR WEATHER SHIP DELTA, 44°N., 41°W. 
collected from weather ships, these graphs show the monthly 
march of mixed layer depth at Atlantic weather ship stations 
COCOA, DELTA, and ECHO. The features of these graphs are fa- 
miliar. As a result of excess absorbed inSolation over emit~ 
ted radiation, the thermocline forms in late spring and early 
summer and is driven progressively deeper by continued mixing 


during the late summer and autumn, until it becomes indistin-= 
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FIGURE 14. MEAN MONTHLY MIXED LAYER DEPTHS FOR WEATHER SHIP ECHO, 35°N., 48°W. 
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FIGURE !5. MEAN MONTHLY MIXED LAYER DEPTHS FOR PAGIFIC WEATHER SHIP NECTAR, 30°N,, |40°W. 


guishable below h50 feet. Minor variations with latitude are 
Meee. Kagurmeas 15, 16; and 7 give the same type of presenta-= 
tion for three Pacific weather ship stations, computed from the 
BT data on file at the Hydrographic Office. In these Pacific 
data it will be noted that, in the early summer, more than one 
value for the mean depth of the thermocline is shown (figure 16). 


This is the transition period when the old, deep winter thermo-~ 
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FIGURE 16 MEAN MONTHLY MIXED LAYER DEPTHS FOR PACIFIC WEATHER SHIP GOLF, 38°N.,158°W 
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FIGURE 17 MEAN MONTHLY MIXED LAYER DEPTHS FOR PACIFIC WEATHER SHIP PAPA, 50°N.,145°W. 


cline is sometimes still in evidence and at the same time the 
new, shallow thermocline is forming, A mean between the two 
would not only be unrepresentative of actual conditions, but 
would give the impression that there is continuity to the curve 
of mean monthly layer depths. Naturally this is not the case; 
one thermocline is distinctly different from the other. 

be Monthly or luner variations = In working with time- 
series BT data, it was noted that long period fluctuations seem= 
ed to exist with periods suggestive of monthly oscillations, and 
superimposed on the seasonal variations, as illustrated in figure 
18, It was considered that, if such fluctuations actually exist, 
they cannot be attributed to turbulent mixins caused by winds or 
sea. Although turbulence will cause an increase in the depth 
of the thermocline, there is no reason to expect the thermoclime 
to prise gradually when the wind or sea subsides. Instead, the 
original thermocline should remain deep after the turbulence 


by wind and sea, a new mixed 
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FIGURE 18. VARIATION OF THE DEPTH OF THE THERMOCLINE WITH TIME AT WEATHER SHIP DELTA, 44°N., 41°W., 1-22 AUGUST 1950. 


layer should form. Since the wave motion observed does not con=- 
form to this description, some other factor must be at worke 

In an attempt to discover if a 29-day lunar variation in 
the thermocline exists, the most complete set of daily BT obser- 
vations available for one year was plotted. These data were for 
the year 19h8 at the Pacific weather ship NECTAR. An attempt 
was made to eliminate the shorter period variations by plotting 
the mean value for éach day rather than each single observation. 
This method was not successful because too many days had only 
one or two observations instead of the expected six; conssquent- 
ly, it was impossible to tell in which phase of the daily oscil- 
lation this single observation was taken, Since short period 
variations still appeared on the curve, a regression curve was 
computed to eliminate them, as shown alin! se alvenbligts) UBS) Gh Biel telah ahs 
graph the seasonal variation already discussed is apparent and, 


in addition, a suggestion of ridges and troughs appears in the 
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curve with periods of approximately one month. Fourier analysis 
was applied to these data in an attempt to discern a true 29-day 
eyeles the results were nescative. This fact. is not surprising, 
since so many other parameters affect the depth of the mixed 
layer that any monthly variation must be treated not as a total 
cause in itself, but rather as a partial cause together with 
other causes. This technique will be discussed later. 

c. Semidiurnal and diurnal tidal variations - It has already 
been mentioned that previous investigators have concluded that the 
greatest short period variations in the depth of the thermocline 
have their origins in an internal wave of tidal character. For 
this reason the main emphasis has been placed on attempting to 
discover techniques which would predict this effect. 

The simplest technique one can use in this case is a purely 
statistical one, such as is used in predicting the tides theme 
selves. It was noticed that some of the time-series traces 
of the thermocline variations at ship stations in the Hastern 
Pacific bore a marked resemblance to the tidal curves at stations 
on the West Coast of the United States. Some of the best data, 
therefore, were given to the U. S. Coast and Geodetic Survey 
with a view to having the data analyzed for the presence of 
semidiurnal and diurnal constituents and eventual ~rediction of 
the depth of the thermocline by means of constants. It was 
fully realized that the depth and intensity of the thermocline 
have an effect on both the amplitude and phase of the wave, but 
PMswooped tiopebae wesulus would Be Wseheuly mor ay pamb euilene 


place’ by USing a digterent set of constants fer each location. 


FIGURE 20. OSCILLATIONS OF THE THERMOCLINE AT WEATHER SHIP DELTA, 44°N., 41°W., AUGUST 1950. 


The analysis of the Coast and Geodetic Survey, however, proved 
negative. It was their opinion that there were deficiencies in 
the data, since, with the number of daily observations varying 
from only one to six, one cannot be sure that the peaks and 
troughs of the waves are clearly defined. There is no doubt 
also that the effect of the depth and the vertical density 
gradient upon the change in phase and amplitude of the wave 
with depth was also a major contributing factor to the failure 
of such an attempt. Illustrations of the type of data which 
were analyzed in this manner are given in figures 20 and el. 
The irregularities in the amplitudes would lead one to believe 
that the point of the Coast and Geodetic Survey is well taken. 
For a complete statistical analysis, observations should be 
taken at much more frequent intervals, in order that the curve 
more closely approximate a continuous trace and all the ridges 
and troughs be delineated. 
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mentioned periodic motion, there are known to exist in the 
oceans internal waves of periods of the order of magnitude of 

a few minutes. The most complete treatment of these short 
period internal waves has been given by Ufford (1947 a, b; ep hes 
These papers derive relationships civing the amplitude, period, 
and wave length of these waves for several special cases of 
multilayered systems. Ufford's work seems to give promise of 
being of practical value in the study of these waves, although 
as yet the present project has not been extended beyond consider= 
ation of the semidiurnal and diurnal periods, Data have been 
collected, however, which illustrate this effect, and examples 
are given in figures 22 and 23. It will be noted in figure 22 
that superimposed on what is possibly the major semidiurnal 
tidal variation are many small oscillations of a period of a 


few minutes. 
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FIGURE 22. SHORT PERIOD FLUCTUATIONS OF THE THERMOCLINE AT 
WEATHER SHIP BRAVO, 56°N.,5I°W, 2000-22002, 13 
SEPTEMBER 1945. 
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FIGURE 23. SHORT PERIOD FLUCTUATIONS OF THE THERMOCLINE AT 
WEATHER SHIP ALFA, 6I°N.,33°W., 1500-1700Z, 22 
SEPTEMBER 1945. 


There is little doubt that short period oscillations as 
described above is one of the main factors that cause discrepan= 
cies between theoretically computed and actually observed sonar 
conditions. However, the limitations of the available data have 
made it impossible to approach this problem in conjunction with 
the longer period waves. It is hoped that in the future arrange-= 


ments can be made for such data to be collected and the problem 


investigated. 


e@. Random variations due to mixing = Since the mechanisms 


that cause mixing in the ocean are mainly meteorological in 
nature, this problem has been discussed in the earlier portion 
of this report. It was concluded that there were two main proc-= 
esses that affect the depth of the thermocline: convectional 
mixing due to instability, and turbulent mixing due to wind and 


sea, When the random effects of these two processes are superim= 
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posed on the periodic variations discussed above, a highly ir- 
regular pattern results. When these data are subjected to an 
empirical treatment, there are two difficulties which must be 
overcome. One is the problem of determining the proper lag 
factor to employ, because the thermal structure does not reach 
equilibrium with the sea and wind instantaneously. The other is 
the fact previously mentioned, namely, that while turbulent mix- 
ing causes the thermocline to get deeper, ceasing of the turbu- 
lence does not necessarily cause the thermocline to become shal- 
low. 

f,. Empirical investigation of thermocline depth - With 
all the effects described above existing at the same time and 
in view of the limited data available, it would be difficult 
to correlate any one of these single parameters with the varia-= 
tion in the thermocline because of the effect of the other fac- 
tors. This procedure has been attempted without success in cor-= 
relating the depth of the thermocline with sea and wind conditions 
alone. 

On the other hand, on the basis of our present knowledge of 
the thermal structure, it can be shown that it would be just as 
hopeless to attempt to include all of the parameters in one re= 
lationship. This would amount, for all practical purposes to 
predicting the depth of the mixed layer for one particular time. 


Such a relationship might be written as follows; 
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In this expression D represents the seasonal effect, and 
may be taken as the mean depth of the thermocline for a particu- 
lar month. The 2nd through 6th terms represent periodic effects 
due to the lunar, diurnal, and semidiurnal waves. In the last 
term, W is some function which represents the mixing processes 
that take place. It might be a complicated function of wind and 
sea, and only experiments will determine its exact nature. ‘The 
contribution to the amplitude of the oscillations of the thermo- 
cline made by the lunar wave would be represented by-j A+B? 
and the phase change rrom the lunar at the thermocline depth is 
given by ten The same relationships obtain for the 
semidiurnal and diurnal waves. The normal procedure in such a 
study would be to evaluate from known data the constants Aj, As, 
By, Bo, etc., and to assume that these constants hold with time 
iomkiat particular place. Assuming this to be’ so, one could, 
therefore, knowing D and the other factors, evaluate D for a 
particular time. 

There is no evidence that either the amplitudes or phases 
of the constituents waves are conservative with time. On the 
contrary, there is evidence to support the fact that both the 
eiuikt senles and phases vary with time and also with depth at one 
particular time. Another difficulty lies in the nature of the 


funetion of the wind which must be used to evaluate turbulent 
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mixing. Ih the light of these facts, it would seem that, based 
on our present knowledge of the present processes involved, a 
precise empirical prediction of the thermocline depth is not now 
possible. 

However, it is possible to make approximations. One of the 
things that we know about the diurnal and semidiurnal oscillations 
is that their amplitudes are dependent to a large extent on) the 
depth at which these oscillations occur. For approximation 
purposes, therefore, the problem may be devided into two parts: 
(a) to predict the mean depth about which the thermocline oscill- 
ates daily and (b) to predict the amplitude of the oscillations 
about this depth. 

The first of these is partly seasonal, partly random, and 
may be partly monthly. The problem involved is essentially one 
of extrapolating the curve shown in figure 19. Attempts have 
been made to devise a successful technique based on such para- 
meters as seme lstieney. surface temperatures, sharpness of the 
thermocline, and a departure of the depth of the thermocline 
from normal. As yet a relationship which can be considered 
sufficiently accurate has not been devised, but work on a success= 
ful method of approximation will be continued in the near future. 

The second problem mentioned above has been attacked with 
some success. First, at several locations in the oceans at 
which many time-series BT observations were available, the 
mean amplitudes of the semidiurnal and diurnal oscillations 
were computed for each of several categories of depth 


at which the thermocline was found. The results are shown 


50 


DISTANCE FROM TROUGH TO GREST (FT) 


ALINIOIA ONY “McOblNoOS ‘YVLOSN dIHS YSHLVSM d1510Vd 


‘SBLVTTINSO 11 HOIHM LNOBV Hid3d NVAW 3H ONY 
S3AVM TWNYSLNI WWOIL SHL SO SGNLIIdWY 3HL NS3SML38 dIHSNOILV13Y @ 


(0) 


S| 


02 


ce 


cr 


Sb 


OS 
6lb-O9E 6CF-O0f 66¢0b2 6£2-08! 62I-Od) 61l-09 6S0 


(13) SNIIOOWYSHL JO H1d30 


“SHLd30 ANITIOOWYSHL ATIVG NVSAW JO NOILNGIYLSIO S9VLNADYSd Wee AYNSIS 


6lv-O9¢ 6GE-00E 662-Ob2 6£2-081 621-02) 61l-09 6S- 
(13) SNMOOWYSHL 40 Hid3d 


08 


sasv9 4O 1N39Y3d 


51 


DEPTH OF THERMOCLINE (FT.) DEPTH OF THERMOGCLINE (FT) 
0-59 60-119 120-179 180-239 0-59 60-119 120-179 


180-239 


PERCENT OF CASES 
ow 
(e) 
ol 
fo} 
(14) 1S349 OL HONOYL WOYS JONVLSIC 


FIGURE 25A. PERCENTAGE DISTRIBUTION OF MEAN B. RELATIONSHIP BETWEEN THE AMPLITUDE OF THE 
DAILY THERMOCLINE DEPTHS. TIDAL INTERNAL WAVE AND Ti1E MEAN DEPTH 
ABOUT WHICH IT OSCILLATES. 


ATLANTIC WEATHER SHIP COCA, 52.5°N., 35.5°W., AND VICINITY 


in figures 2) and 25. In each figure the graph at the left indi- 
cates the percentage of the total observations which showed the 
mean depth of the thermocline at each of the depth categories and 
the grapk at the right indicates the mean amplitudes plotted 
against each of the depth catesories. It is evident that there is 
a direct relationship between the mean depth about which the 
thermocline is oscillating and the amplitude of the oscillations 
themselves, It will be noted that variations from a straight 
line relationship exist. In particular it was noticed that the 
slight decrease in the amplitude at the sreater depths shown in 
the right-hand graph of figure 25 was also characteristic of 

some of the data computed for other areas of the oceam but. not 
presented in this report. It is not known at present whether 
this is due to the relatively few observations taken where the 


thermocline was found at the Sreater depths or whether there is 
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some physical reason as yet uhlmown, In general, however, all 
areas in the ocean showed an increase in amplitude with an in-=- 
crease of the mean depth of the thermocline, which agrees sub- 
stantially with the two-layer theory. 

Such depth versus amplitude data may be used to devise an 
empirical relationship which will be applicable to one particu- 
lar place in the ocean. As mentioned before, appendix A pre- 
sents a derivation of an expression for the amplitude of a wave 
along the interface between two fluids using the linear two 
layer theory and taking into account both the rotation and the 


spherical nature of the earth. This amplitude was found to be 


Ko Poe K9PrFIA Kg Py | 
_ Reos¢ h, RrpeF dARopyF «Rpg Cos¢ 


 RRcos¢ ARcosé/KgeFA, egerF KaPy 
nN hy ARop,F Rrp,F Regcos¢ 


(20) 


The nature of the terms is explained in appendix A. It will 

be seen that many of the terms are functions of the latitude, 
radius of the earth, and the wave length. If we may assume, 

as Haurwitz did, a constant wave length for the semidiurnal 

and diurnal waves, then all of these terms become constant for 

a particular place. In addition, for deep water it is not 
possible for h,, to change very much relative to h;, since it is 
so much greater in magnitude. Since p, for all practical purposes 
may be considered constant, it is then possible to include both 


the depth and density of the lower layer in the constant terms, 
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By grouping the constants for each term together, the expression 


then becomes 
—— 
syn Py 2) 


Here the Ks are arbitrary constants to be determined, Re- 


arranging terms and grouping the constants once again, we have 
OG CarerGs + Gaher Capet =i (22) 


The expression is now in a form which may be subjected to 

a curve fitting procedure. From analyzed BT data for a parti- 
cular area, such as is presented in figures 15 and 16, the 
amplitude and the depth of the thermocline hy may be determined, 
and a good approximation of p. may be arrived at from the sur-= 
face conditions. By using sufficient data, the constants in 
equation (13) may be computed. This was done for the data 
available in the vicinity of the Pacific weather ship NECTAR 
(the data from this area ane summarized in figure 2). When the 
constants were evaluated and the resulting expression solved for 


the amplitude, the following empirical relationship was obtained 


_1-97.443 h,+ 95.0908 p, h, 


= 2 
241= 3800824 +273.3448 P, cr 


It was more convenient in processing these data to use the 
distance from trough to crest, rather than one-half distance. 


Therefore, in this expression 26, represents not the actual 


54 


amplitude but twice that value. 

This equation for the present should be considered a rough 
approximation only, since the actual data were so scarce and 
so infrequent that a true approximation of a continuous curve 
was not very often obtained. In order to test the equation, a 
few observations were available from weather ship NECTAR other 
than those used in computing the constants, A listing of the 


computed and observed heights (twice amplitude) follows: 


Computed Observed 
28 26 
el 15 
S) 36 
Be 27 
105 8 
vB) 2 


The correlation coefficient between the above computed and observed 
heights is 0.69, which, for only six cases, is of no statistical 
significance. The results, however, seem to indicate a reasonable 
basis for further work on this type of approach. The information 
contained in equation (23) may be summarized by a single nomogram, 
as shown in figure 26, For a given density of the upper 
layer and for a given mean depth of the thermocline, the am- 
plitude of the wave along the thermocline may be read off. 
An unusual feature of this equation is the fact that there 
exists a thermocline depth where the amplitude is independent 
of density, and remains nearly constant. No particular physical 
interpretation of this feature can be given. It should be 


emphasized once again that equation (23) and figure 26 are 
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based on data taken from one location; therefore, they are 
Molt Only in that area. 
3. ‘The continuous density model 

So far the investigation has centered around a model 
in which the wave occurs along the interface between two homo- 
geneous layers of water. If more than one interface are present, 
then several internal waves can result. Concerning this event- 
uality, the Oceans (1932) says "...the greater the number of 
boundary surfaces the greater the number of internal waves. On 
the basis of this reasoning, when the density varies contin- 
uously with depth one should expect an unlimited number of 
possible internal waves." The effect of this concept upon an 
investigation of periodic variations in the oceans frequently 
would be to give results different in nature from those fiven 
by the two-layer system. With many infinitesimal vertical 
changes in density occurring simultaneously, it does not 
necessarily follow that the largest amplitude will be found at 
the thermocline. A complete statistical analysis of this 
problem has been presented by Seiwell (1942). This consists of 
a complete harmonic analysis of the oscillations at several depths 
at "Atlantis" anchor station 3245. One of the main conclusions 
of this work was that the smaller the vertical density gradients 
at a given depth, the larger will be the amplitudes of the 
oscillations, and the smaller will be the temperature chanses 
with time at that depth. In the case of the particular data 
that were analyzed in this study, it was found that the greatest 


amplitudes occurred at a depth somewhat greater than that of 
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the thermocline. 


The picture thus presented by the continuous density model 
is one much more complicated in nature than that of the two-layer 
system, We are now concerned with an ocean in which the ampli- 


tudes and phases of the internal wave are continuous functions 
of depth, 

The most complete theoretical treatment of this complicated 
problem has been given by Fjeldstad (1933). Assuming free waves 
in a medium in which density varies continuously with depth and 
neglecting the effect of friction, Fjeldstad derives the follow- 


ing second order differential equation 


dw 4 =O 
dz? g owe (2h) 


with boundary conditions at the bottom (w#0, Z=0) 

and at the surface (w=0, Z=sh), 

In this equation, w represents the relative displacement of a 

water particle from its equilibrium position. This means that 
these amplitudes are not unique, but only that they are deter-= 
mined in such a way that the ratio of true amplitudes between 

any two fixed depths is equal to the ratio of the w's for the 

same depth. The term ¢is taken as the internal stability of 


the water column, or 


| dp 
p=— SP 2 
5 ab (25) 
2 ° 
hig 1s an unknown parameter which depends on the vertical 


density distribution. 
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it is necessary in general to carry out this integration 
by numerical methods, and an infinite number of solutions corre- 
sponding to an infinite number of internal waves are possible, 
This integration yields the solution for the relative amplitudes 


Wa 


The internal wave motion is then assumed to be of the 


following form 


n=A cosct+Bsinct (26) 


Where the computed aniplitudes may be evaluated by the following 


summations 


A= 3a, B= b, Wy (27 a,b) 


The coefficients ay, and b, may be evaluated by the use of 


the Fourier theorem 


h h 
ane k [Ad 2% dz Dek |B dc2rm dz (298 a,b) 
{0} {e) 


It will be noted that in computing the Fourier coefficients a, 
and b,, it is necessary to have observed amplitudes A and B in 
equations (28). These coefficients are then used to compute the 
actual amplitudes in equations (27). Wo matter what the nature 
of the orthogonal function w, the Fourier theorem requires that 
the observed and computed amplitudes must be identical, if an 
infinite number of orders are used, However, Fjeldstad demon- 


strates that, by employing his particular W, good agreement is 
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obtained between the observed and computed amplitude in four 
orders. This is a very strong mathematical indication of the 
validity of his theoretical formulation. His differential samae 
tion, therefore, has been widely employed in studies by other 
authors, notably Seiwell] (192) and Lek (1938). 

It must be emphasized that Fjeldstad was concerned prima- 
rily with validating his theory, and not with computing ampli- 
tudes at a time and in a place where no amplitudes were known. 
His work has been examined in the present investigation with a 
view to projecting his results into the future in order to pre= 
dict the character of the semidiurnal and diurnal internal wave 
in a given place at some future time. It was hoped that in his 
concepts could be found some conservative or predictable factor 
which could be employed in a technique for successful predic- 
tione 

The most convenient and appealing idea was that the Fourier 
coefficients themselves were conservative with time for any one 
place. However, no satisfactory theoretical defense for this 
thesis could be found, since it seems clear from the above equa-= 
tions that, if the density distribution changes with time, not 
only would the amplitude A change but also the coefficients an. 
This does not mean, however, that this concept as a possible 
basis for a prediction technique should be discarded entirely. 
In cases of slow time rate of change of density distribution, it 
is possible that the coerficients would remain reasonably con-= 
stant with time; from a practical standpoint this would be almost 


as useful as absolute constancy. It may be that even for large 
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density changes with time there would not be much change in the 
coefficients, but this has not as yet been shown. However, 
Pekeris (1937) and Wilkes (1949) have investigated a similar 
problem in the atmosphere, with the difference that they postu- 
late a driving force rather than assume free waves. An exami- 
nation of their work reveals the hope that a similar treatment 
of the oceans will show that the Fourier coefficients will re- 
main relatively conservative with time, since analogous terms 
in the work of Pekeris depend only on the vertical density dis- 
tribution and the location. 

Since the approximate conservativeness of the Fourier co- 
efficients could not be demonstrated, it was decided to proceed 
upon the assumption that a change in these values occurs. It 
was further assumed that when this change occurred, it occurred 
in such a manner that the change in she time components of the 
semidiurnal and diurnal waves was a minimum. This assumption 
was predicated on the concept that, for any new set of initial 
conditions, the ocean like all fluids will reach equilibrium 
with these new conditions with the least possible change in the 
motion. 

Many possibilities for minimizing changes other than the 
above present themselves. For instance, one could attack the 
problem of minimizing the change in energy from one set of 
conditions to another by minimizing the change in the square 
of the amplitude. Another possibility would be to attempt to 
minimize the change in the amplitudes and the phases themselves. 


While these technicues obviously have merit, they were super- 
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ficially investigated and found to require a computational 
treatment far longer than the limited staff of the project 

could supply. With the aid of machine computation, however, 
this problem could be attacked in the future. Therefore, it was 


decided to proceed with the original assumption. 


Let Aj be one of the two constituents of the amplitude of a 
wave computed from the observed data at one particular time ac- 
cording to Fjeldstad's theory. Let A> be the corresponding un-= 
known component for some future time at the same place. At this 
future time, let one observation of vertical density distribution 


be available so that the w's are known, 
4 2 
Let = (Gin Win ~ConWon) =I (don) 
Misia 2 (29) 


Now by minimizing the function I we are minimizing the change of 
the component of the amplitude of the wave Aj (see equation 27). 
The minimizing of the function I may be accomplished by the so= 
lution of four simultaneous equations which will yield the un-= 
net Aone This computation is long and involved and, for the 
example mentioned below, will not be repeated here. It was car= 
ried out by desk calculator by using the inverse matrix approach. 
For the harmonic analysis and the numerical integration, tech- 
niques identical with those of Fjeldstad were employed. 

The type of data required to test this method is that ob- 


tained at a standard anchor station of at least 2: hours! dura- 
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tion, repeated ae some future time in exactly the same place. 
Adequate data such as these are not as yet available. The only 
repeated anchor station the investigators could find in all lit- 
erature was that taken in the Florida Straits by two U. S. Navy 
Hydrographic Office Survey ships, the SAN PABLO and the REHOBOTH, 
on the th and 9th of December 1950. The deficiencies of these | 
data were many, and they will be discussed later. However, since 
these were the only data available, it was decided to employ them, 
more as an instrument for developing a technique than es an hon-= 
est test of the method. 

The data for h December were analyzed according to Fjeld- 
stad's method and the values for the Fourier coefficients are 
presented in table VIII. By numerical integration the relative 
amplitudes, w, were computed for the lh December data at 5 chosen 
depths. These values are given in table IX. The data for 9 
December were analyzed in the same manner and the relative amp= 


litudes are presented in table X. 


Table VIII 
Order Diurnal Semidiurnal 
- .980 oO 
eis tae a 
2h - . 6.8 22930 


Fourier coefficients as computed from 
anchor station of | December 1950. 
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Depth 
(meters) 


Relative amplitudes, w, as computed from 
anchor station of lh December 1950. 


Depth 
(meters) 


Relative amplitudes, w, as computed from single 
density distribution on 9 December 1950. 
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Table XI 


Semidiurnal 


Predicted Fourier coefficients 
for 9 December 1950 


Table XII 


Depth Ay 
(meters) Predicted 


Predicted and observed absolute amplitudes 
(meters) of the diurnal internal wave on 
9 December 1950 
The data presented in tables VIII, IX, and X may now be 
serted in equations (29) and the solutions for the predicted 


Fourier coefficients on the 9th of June obtained. These are 


eealyeim alia iallolle Sir By using these predicted Fourier coeffi- 
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in- 


cients, the data for the relative amplitudes for the 9th of De- 
cember given in table X, and equation (27a), it is now possible 
to predict one of the components of the internal wave for the 
9th of December. The values of this component for each of the 
chosen depths are given in table XII. Finally, as a verifica- 
tion, it was necessary to compute by harmonic analysis the ac- 
tual observed components of the amplitudes. These observed am- 
plitudes are also shown in table XII. In like manner it is pos- 
sible to compute the other component of the amplitude, 

A comparison of the predicted and observed amplitudes as 
given in table XII indicates that the results, while promising, 
are not entirely satisfactory. In general, it is considered 
that there are two main reasons for the discrepancies in the 
forecast. First, there are deficiencies in the data themselves. 
The anchor stations on the lth and 9th of December were each 
taken by two ships anchored approximately 30 miles apart. lHach 
ship took 6 observations during the 2i-hour period of the th 
and these were combined to form one set of l2 observations for 
that one day. The procedure was repeated on the 9th. The fact 
that the ships were anchored so far apart undoubtedly intro- 
duced reasonably large errors in the harmonic analyses. Further= 
more, the anchor stations were located in shallow water and in 
an area in which reasonably strong currents exist. All these 
factors tended to contaminate the results. Second, there is a 


possibility that the incorrect parameter was minimized. ‘his 


has already been discussed. lia wae uliouies Wwe 1s jolleinmsc! we 


conduct other trials minimizing other parameters. 
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Another interestins feature of these results is the fact 
that the prediction was much more accurate at greater depths 
than at the surface. The example shown in table XII includes 
only one component of the diurnal wave; actually both components 
for both the semidiurnal and diurnal wave were computed. In 
all four cases the results displayed the same inaccuracy near 
the surface and the same accuracy at the greater depths. Whether 
this particular inaccuracy is due either to one of the above 
two reasons or to some near-surface phenomenon not considered 
in the computations has not as yet been determined. 

It will be noted that the technique outlined above does 
not, according to strict terminology, constitute a true fore- 
cast, since either an observation or an assumption of a present 
density distribution is necessary to employ the method. However, 
Dee's Penoricened a large step forward to be able to estimate 
the amplitude of the wave from a single observation, rather than 
from a complete enietiote beaten with many time-series observations. 
If some such method as described above were to be developed, it 
would automatically increase the chances of predicting the 
internal wave, for it is not possible to occupy a complete anchor 
station every time a forecast is desired. If one harmonically 
analyzed all existing and future anchor stations and computed 
the Fourier coefficients for each, then at any future time the 
internal wave could be estimated for each of these areas. The 
only data one will need for this operation would be a single 
density distribution. This observation could be either (a) @ 


recent observation talken in the area, (b) an estimated distri- 
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bution based on the latest available data and extrapolation, or 
(c) a predicted distribution based on climatological data for the 
particular month desired and tempered by the latest information 
and by experience. It is believed that, if a mass of data as 
mentioned above was collected and such a technique developed, 
it would be of great value in delineating one of the main effects 
upon the thermal structure of the ocean. With this effect evalu- 
ated, it would then be possible to proceed to the study of other 
factors. 
Iii. Data 

The single feature that has had the greatest effect on the 
course of this investigation has been the lack of data suitable 
for conducting it. When the approach to be employed in this in- 
vestigation was first decided upon, the type of data that would 
be necessary to employ this approach immediately became clear. 
It was evident that most of the data that were available from 
various sources had two serious defects: (a) The data wereobtain- 
ed by the classical technique; that is, both time and distance 
were allowed to vary in obtaining the serial observations, and 
(b) none of the data contained measurements of a sufficient num-= 
ber of parameters to evaluate completely all the effects on the 
thermal structure. 

The first of these defects is by far the more important, 
It is the opinion of the investigators that, if problems such as 
this are to be solved eventually, more attention must be given 


to eliminating the factor either of time or distance from the 


collection of data, One of the first things that becomes evi- 
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dent in recent oceanographic literature is that, when analyses 
are made either of time-series or synoptic data, the ocean is 
found to be much less invariant with time than was heretofore 
assumed. No better illustration of this thesis could be given 
than that of Defant (1950) who demonstrated that, when one elim- 
inates the effect of each observation being taken at a time and 
place different from every other observation, an entirely dif- 
ferent and more credible pattern of the oceanographic variables 
results. In this case Defant eliminated the effect of the tidal 
internal wave, but undoubtedly a thorough investigation of other 
factors, such as the effect of horizontal gradients in the ocean, 
would give the same result. 

It was clear that what was needed for this investigation 
was a series of repeated anchor stations taken at one or more 
HocatLtons. | these anchor stagions initadally should be in 
locations where the variables which could not be accurately 
measured would have little or no effect, such as an area where 
relatively little current exists and where there are no abrupt 
changes in bottom topography. These observations should be 
repeated every hour or two, and frequent BT's taken at regular 
intervals. They should also include simultaneous observations 
of such other parameters as might be important to changes in 
the thermal structure, such as incoming and outgoing radiation, 
transparency, air temperature and humidity, wind, and sea and 
Swell. It.was considered that not too much could be learned 
from single anchor stations of one day's duration since they 


afford little or no opportunity for testing techniques or 
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conclusions. 

The alternative to anchor station type of data would be 
a series of observations taken simultaneously by several ships 
each at a different location. This type of synoptic data would 
be very useful if one wished to consider the effects of hori- 
zontal advection on the thermal structure, since it would give 
an estimation of horizontal pradients. Actually the optimum 
would be a combination of the two types of data. Except for 
the possibility of expanding the weather ship program, however, 
the cost of any such network as a regular procedure would be 
prohibitive. 

This then was the type of data that was desired. Searches 
were made in all possible sources in order to obtain the data 
which best fit the requirements, but with very little success. 
By far, the bulk of the time-series data available consists of 
BT data taken from weather ships. In addition the only repeated 
anchor stations that could be found were the set of observations 
used in the example in the text. The deficiencies of boven) eiesm 
types of data have been discussed above in connection with the 
discussion of the use to which the data was put. 

Although the data acquired in time to be of use in this 
report were far from adequate, much progress has been made both 
in arranging for complete collections of suitable data and in 
the development of new observational techniques aimed at increas- 
ing the value of the data, Much of this work has been carried 
on in the Hydrographic Office itself. Arrangements have been 


made for future cruises of the Hydrographie Office survey Ssiiaps 
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to include among their observations long term anchor stations in 
suitable locations with closely spaced serial observations of 
all measurable parameters. One of these cruises, AMOS XI, is 
now in progress and the data will be available in the very near 
future. In addition, in cooperation with the Office of Naval 
‘Research, plans have been completed for conducting a series of 
cruises of five days! duration in each month to collect similar 
data to be utilized synoptically with the observations taken at 
weather ship HOTEL. The U. S. Navy Hydrographic Office has been 
very active in improving the instrumentation required for this 
type of oceanographic data. The radiation data presented in the 
first section of this report constitutes one of the first sets 
of data from a recently developed pyrheliometer array now in- 
stalled on the U. S. Navy Hydrographic Office Survey ship REHO= 
BOTH. New photometric equipment is now being tested with a view 
to improving methods of obtaining transparency data with depth. 
In addition, the Hydrographic Office has developed to a stage of 
near completion a thermopile to supplement BT observations which 
is expected to yield much valuable information on the microther= 
mal structure of the ocean. This device is of great interest to 
the present investigators, not only because of its obvious value 
in estimating very small changes in the thermal structure due to 
processes such as radiation and evaporation, but also because of 
its great possibilities as an instrument for measuring fluctua- 
tions in temperature of a periodic nature. Such an instrument, 
when held at a constant depth, would give a good indication of 


these periodic fluctuations; several such units, each at a dif~ 
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ferent depth, would delineate the relationship between amplitude 
and phases at these different depths. 
I¥. Recommendations for future work 

The project as it has developed so far has opened several 
avenues of approach which, if followed, should yield interesting 
and useful results of a practical nature. These consist partly 
of extending approaches already employed by utilizing more ade- 
quate data, and partly new approaches suggested by various phases 
of the work and by the literature. They may be summarized as 
follows: 

ae A re-evaluation shovld be made of the first portion of 
this study dealing with the direct effect of meteorological pa- 
rameters, such as radiation, evaporation, and wind. It seems 
clear that the small changes in the thermal structure caused by 
changes in these parameters must have an important effect on 
operational problems connected with sound. It is believed that 
with increased accuracy of the instruments discussed above, es-=- 
pecially those such as the thermopile dealing with the measure-= 
ment of temperature, much useful information may be obtained. 

b. It is planned to continue and expand the effort to ar- 
rive at useful empirical relationships aimed at predicting the 
mixed layer depth. 

ce Further computations using the continuous density model 
of Fjeldstad are proposed, by utilizing more accurate data and 
making such revisions of the technique as necessary. 

d. It is also planned to investigate the possibility of 


attacking the problem of tidal oscillations in a direct manner, 
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i.e., by considering the tidal force itself. This would be 
a purely deductive method stemming directly from the differ- 
ential equation. 

e. Another possible approach now being considered is an 
2S) amo to develop a purely noncausal technique of predicting 
the thermal structure by using the principles of statistical 
mechanics. ‘This would probably employ the largest possible 
collection of BT time-series data and would be developed by 
using standard techniques of time-series prediction. In such a 
study it would be necessary to assume some particular statisti- 
eal form of the distribution of the components of the periodic 
oscillations. An interesting treatment of this type was applied 
by Pierson (1952) who demonstrated prediction in time and space 
of a wave train whose displacements from the mean displayed a 
Gaussian distribution. 

f. All the BT data at the Hydrographic Office are now being 
analyzed for the occurrence of synoptic situations, or situations 
in which there was a reasonably large horizontal distribution 
of observations for a given period of time. If good examples of 
such distributions are found, it is intended to attempt to apply 
well-known techniques of extrapolation and kinematic analysis, 


with a view to developing a useful prediction tool. 
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